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New Neighbours: Modelling the Growing 
Population of 7-ray Millisecond Pulsars 

C. Venter, A.K. Harding, and TJ. Johnson 



Abstract The Fermi Large Area Telescope, in collaboration with several groups 
from the radio community, have had marvellous success at uncovering new 7-ray 
millisecond pulsars (MSPs). In fact, MSPs now make up a sizable fraction of the 
total number of known 7-ray pulsars. The MSP population is characterized by a va- 
(~| [ riety of pulse profile shapes, peak separations, and radio-to- 7 phase lags, with some 

Q^' members exhibiting nearly phase-aligned radio and 7-ray light curves (LCs). The 

MSPs' short spin periods underline the importance of including special relativistic 
$-H ■ effects in LC calculations, even for emission originating from near the stellar sur- 

c/3 I face. We present results on modelling and classification of MSP LCs using standard 

^, ■ pulsar model geometries. 



1 Introduction 



To date, Fermi Large Area Telescope (LAT) has uncovered more than a dozen new 
y-ray millisecond pulsars (MSPs) 1171 . These MSPs were discovered by folding the 
(^ • y-ray data with the known radio period, using timing solutions provided by various 

t^ I radio telescopes. Discovery of the first eight Ferm/-LAT MSPs 1 1 1 indicated that 

these MSP pulse profiles and spectral properties mimic those of the younger pulsars, 
and that there are non-zero lags between the radio and 7-ray peaks. However, the 
. , discovery of PSR J0034— 0534 ^ exhibiting phase-aligned 7-ray and radio peaks 

rS ' provided evidence for co-location of the radio and 7-ray emission regions. This 

j3 ■ phenomena has previously only been observed for the Crab pulsar Recent detection 

of 7-ray pulsations aligned with the radio peaks 1231 141 from the famous 'first' MSP 
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PSR J1939+2134 0, as well as the first 'black-widow' MSP PSR J1959+2048 IfTSJI . 
now provide two more examples of this exceptional behaviour. Publication of a 
number of new 7-ray MSP detections are awaited 1221 (some resulting from follow- 
up radio observations of unpulsed Fermi-LAiT sources which appear in the First 
Fermi-LP^ source catalogue), and the growing MSP population now constitute a 
significant fraction of all Fermi-LPsI pulsar detections to date 12]. 

The y-ray MSPs exhibit a variety of light curve (LC) shapes, and have been mod- 
elled using standard two-pole caustic (TPC), outer gap (OG), and pair-starved polar 
cap (PSPC) models 1261 , an annular gap model ITTI , lately also altitude-limited TPC 
/ OG models 125JI , and lastly using a force-free magnetic field Q. We elaborate on 
some of these models in Section|2] (See ITTl for a review of y-ray pulsar population 
modelling.) 



2 Pulsar Models 

Below is a brief description of some modelling aspects (see l26l for details). 

B-field As before, we assume the retarded vacuum dipole field ITOl [T3JI to be the 
basic geometric structure of the MSP magnetosphere. The magnetic moment is in- 
clined by an angle a with respect to the spin axis, while ^ is the observer angle. 

Caustics Due to the rapid increase of the co-rotation velocity with altitude, ap- 
proaching c as one reaches the light cylinder, it is important to include relativistic 
effects when calculating LCs. The effects of (i) B-field line curvature, (ii) time-of- 
flight delays, and (iii) aberration combine to form intense caustic emission on the 
trailing field lines, while the emission is spread out on the leading field lines l20l[T4ll . 

Polar Cap (PC) Model The PC model l9| assumes extraction of primary charges 
from the stellar surface, acceleration of these charges to relativistic velocities along 
B-field lines, and subsequent 7-ray curvature radiation and pair creation at a very 
low altitude. A pair-formation front (PFF) is formed which screens the accelerating 
E-field (£11) parallel to the B-field. 



Faded PC Model This is a geometric approximation H251 . where we use a rising 
and falling exponential function to modulate the assumed emissivity close to the 
stellar surface, as motivated by pair cascade modelling. 

PSPC Model When E\\ is not large enough to accelerate primary electrons to suf- 
ficiently high energies so pair cascades will occur, it will be unscreened, and one 
expects acceleration and radiation from the full open field line region above the 

PCQUlIll. 

SG / TPC Model When enforcing the boundary condition that E\\ vanishes on the 
last open B-field lines (LOFLs), the PFF altitude becomes dependent on magnetic 
co-latitude, and narrow gaps form just inside of the LOFLs |5 1. The SG model ETI 
may provide the framework for the geometric TPC model \VI\ . 
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Fig. 1 Altitude-limited TPC / OG LC fits for PSR J1939+2134 (panel [a]: y-rays\ panel [b]: 



radio) with (aTPC^^TPC) ^ (65°, 80°), (aO«,C°°) = (70°, 90°), (R, 



r,y,TPC p/.TPC 



(0.12,1.05)Rlc. (C-R^bSP) 



(0.12,1.00)Rlc, (<,r"^-«S^x°"^) 



(0.93,1.20)Rlc, 

^^radio.OG^^radio.OG^ = (0.75, 1. 10)Rlc- Fits for PSR J1959+2048 (panel [c]: 7- 
rays; panel [d]: rarfio) with (a^P^ ,jtpC) ^ (60°, 80°), (a°°,CO°) = (55°, 85°), 

7,TPC pT.TPCx _ i-pr.OG p^OGn _ /n r t r orvM> / ciadio.TPC nradio.TPC^ 



«;:^R£a^) = «r.«^^) = {o.i2,i.2o)rlc, «? 

'"'"°'°°,RSt'°°) = (0.90,1.10)RLC. 



,R|r 



= (0.95,1.05)Rlc, 
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min 



OG Model The Goldreich- Julian charge density lfT6ll changes sign at the null charge 
surface (NCS), and OGs are assumed to form close to the LOFLs, above the NCS [S] 

ED. 

Altitude-limited TPC / OG Models We lastly investigated geometric TPC / OG 
models with arbitrary lower, ^min. and upper, ^max. radial gap extensions ||3|- 



3 Results 



First 8 MSPs We found |126| that 6 of the 8 Fenni-LAT MSPs may be modelled 
using standard TPC / OG models, while 2 are exclusively fit by PSPC models. 

PSR J0034~0534 We found altitude-limited TPC / OG LC fits 0, with (a, Q = 
(30°,70°), (/?L,/?Lx) = (0.12,0.9)/?LC, and (/?;^f^°,Cfj°) = (0.6,0.8)/?lc, with 
Ri^c = cP/2n, and P the pulsar period. There also exist solutions ll25l for the faded 
PC model for (ci-X) = (10°, ^ 34°), although the first ones seem more credible. 
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PSR J1939+2134 and PSR J1959+2048 Figure [1] indicates altitude-limited TPC / 
OG LC fits with aligned radio and 7-ray pulses ||4|. 



4 Discussion and Conclusions 

Our results indicate that the high-energy (HE) radiation comes from narrow gaps in 
the outer magnetosphere, requiring the presence of robust pair formation to screen 
much of the open field line region. However, MSPs do not produce enough pairs in 
the standard vacuum dipole models to screen Eu. Standard geometrical models still 
provide suitable LC fits though, and the resulting a and ^ are in reasonable agree- 
ment with fits from radio polarization measurements. In the case of phase-aligned 
7-ray and radio pulses, we infer the existence of co-located emission regions, as well 
as high-altitude radio emission. Future work includes investigations of HE spectra 
using full radiation codes, collective emission from MSPs in globular clusters, and 
pulsar population studies. 
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